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Abstract 

We study the most promising signals of Constrained Minimal Supersymmetry de- 
tectable at a luminosity upgraded 2 TeV Fermilab Tevatron collider. Using a full event- 
level Monte Carlo based on Pythia/ Jetset, we simulate the trilepton signal examining 
in detail the effect of constraints on the parameter space. We also simulate the monolep- 
ton and dilcpton signals, the $t + jets signal, and the signals of stop production in 
supersymmetry all with full Standard Model backgrounds with realistic detector cuts. 
We find that large fractions of parameter space can be probed (or eliminated if no 
signal is found), but mass limits on charginos and neutralinos are not possible based 
solely on the trilepton signal. Detection efficiencies depend strongly on supersymmetry 
parameters beyond simply the neutralino and chargino masses; analyses (experimental 
or theoretical) that do not include this will draw misleading conclusions. Finally, we 
comment on how searches at LEP II will complement searches at Fermilab. 
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1 Introduction 



The search for supersymmetry (SUSY) should be one of the primary goals at the Fer- 
milab Tevatron and any future collider planned or in construction. While e + e~ colliders 
generally have the advantage of unambiguously finding or excluding superpartners up to 
m ~ their mass reach is strictly energy limited. LEP II, the highest energy e + e~ col- 
lider approved, will be able to probe chargino masses up to about m\y but no further. A 
useful complement to LEP II's important contributions to the search for supersymmetry 
is provided by a luminosity upgraded Fermilab Tevatron pp collider at y/s = 2 TeV. As 
we shall subsequently describe, if the correct supersymmetric theory has a chargino mass 
below mw, the upgraded Tevatron might possibly fail to discover it (hence the importance 
of LEP II); however, a high luminosity Tevatron does have the capability of discovering 
supersymmetry at mass scales far exceeding the capabilities of LEP II. 

In this paper we report the results of our simulations of many supersymmetric signals for 
a pp collider at y/s = 2 TeV. To this end, we introduce an implementation of all tree level 
MSSM processes and decay modes in the event generators Pythia and Jetset |2|, [|| . This 
implementation of the MSSM is not only a cross check of ISASUSY 0] (see also ||), which 
uses the ISAJET || generator, but includes the refinements inherent to the Pythia/ Jetset 
system. The full details of the event generator will be published in a forthcoming technical 
report [fjj. 

An equally important task is to run the simulations on supersymmetric solutions which 
could be the correct theory of nature. Therefore, we have simulated supersymmetric events 
in the context of the Constrained Minimal Supersymmetric Standard Model (CMSSM) ||. 
Briefly, the CMSSM enforces gauge coupling unification at the "unification" scale ~10 16 GeV, 
enforces proper electroweak symmetry breaking, assumes common scalar masses (mo), com- 
mon gaugino masses (mi/2), and common trilinear scalar soft couplings (Aq) at the unifica- 
tion scale, requires R-parity conservation, and imposes all known experimental constraints 
such as limits on b — > sj decays, invisible width constraints onZ^ & t LEP, etc. Each 
supersymmetric "solution" , with all of its masses and mixings, is determined by five input 
parameters 

m ,m 1 / 2 , A), tan/?, and sign(^,). (1) 
(We use nit = 170 GeV.) Our convention is that tan/3 = v u /vd, where v u (vd) is the vev 
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which gives mass to the up-type (down-type) fermions. Our convention for the relative 
sign of ji is reflected by our choice of the chargino mass matrix element X22 = — M an d the 
neutralino mass matrix elements Y34 = Y43 = +/i (which is the opposite convention from 
Haber and Kane ||]). 

The requirements we impose on our supersymmetric solutions are consistent and in- 
terrelated. While it is not absolutely necessary for nature to follow all of our theoretical 
assumptions (such as common scalar masses), we do wish to emphasize that each require- 
ment supports the other requirements to some degree. For example, the experimental 
requirements such as limits on flavor changing neutral currents support, although do not 
absolutely require, the theoretical preference in supergravity that all scalars have a common 
mass at the high scale. Similarly, gauge coupling unification supports the notion of R-parity 



conservation [10, 11], implying an absolutely stable lightest supersymmetric particle (LSP), 
which is the lightest neutralino Xi- Proper electroweak symmetry breaking, which is metic- 
ulously enforced in the CMSSM, gives the LSP the right properties to be a natural weakly 
interacting cold dark matter particle, explaining many astrophysical observations ^L% ]. 

One of the strengths of the CMSSM is that it is strongly constrained and theoretically 
restrictive. We believe it is remarkable progress that it is possible to construct consistent 
supersymmetric models that incorporate all the above constraints and provide rich phe- 
nomenological predictions. In general it is always appropriate to study the simplest (hence 
minimal) theory that is consistent with what is known from the Standard Model (SM), 
while simultaneously extending our understanding of nature. The CMSSM is minimal in 
the sense that it has the spectrum and group structure of the SM plus superpartners while 
implementing the minimal supergravity boundary conditions at the unification scale. 

Using the CMSSM framework we have performed event-level analyses of the chargino- 
neutralino trilepton signal, the chargino-chargino and slepton-slepton dilepton signals, the 
slepton-sneutrino and chargino-LSP monolepton signals, the squark/gaugino $t + jets sig- 
nal and the signals from stop production and top decays to stop assuming various integrated 
luminosity scenarios (200 pb _1 , 2 fb _1 and 25 fb^ 1 ) for an upgraded Fermilab Tevatron col- 



lider. There have been several previous studies of supersymmetric signals [13, 14] all of 
which are useful and demonstrate the possibility of detecting supersymmetric signals at 
Fermilab, but often make unrealistic assumptions about supersymmetry or make simplifi- 
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cations in the simulation of the signal and background. This study is more general; we find 
several results different from other studies which stem from two basic principles unique to 
our approach. First, we only consider values of SUSY parameters which are consistent with 
the CMSSM. Second, within this general framework, we examine all parameters rather 
than a set of special cases. 

It is possible to construct supersymmetric frameworks that modify some of the theoret- 
ical assumptions and this may be necessary in the future. However, since we sample the 
supersymmetric parameter space with thousands of solutions, it is likely that many alter- 
natives will lie in regions we have already covered in this study. Hence, the conclusions of 
this paper are not likely to be significantly modified. In fact, there are already some hints 
from LEP data that suggest going beyond the minimality of the CMSSM is necessary due 
to the Rb and a s measurements |jlq , |l6| , |l7f . If these measurements are manifestations of 
non-minimal supersymmetry, then we expect very light charginos and stops to exist, and 
therefore the Fermilab Tevatron will have a much greater chance to discover them. 

We organize this paper as follows. In Section 2 we briefly explain the event simulation 
needed for every allowed solution in the CMSSM parameter space. Section 3 contains the 
motivation and results for each signal we have studied. We perform a full background analy- 
sis on every signal, specifying the cuts needed to reduce backgrounds. From the background 
estimates we determine the detectability or significance of each signal for particular inte- 
grated luminosities. In Section 3.1 we examine monolepton and dilepton signals, in Section 
3.2 we examine the trilepton signal and in Section 3.3 we examine in detail the effect of 
constraints on supersymmetric parameter space with trilepton detection as an example. In 
Section 3.4 we examine the I^t + jets signal, in Section 3.5 we examine the signals from stop 
production and in Section 3.6 we comment on the signal from top decays to stop. In Section 
4 we discuss going beyond the CMSSM, including possible effects on supersymmetric signals 
and detection. Finally, in Section 5 we conclude with a summary of our results and a brief 
discussion of how a high luminosity Fermilab Tevatron would complement LEP II. 

2 Monte Carlo Simulation 

Each solution defined by the supersymmetric parameters mo, J7H/2, A)j tan/3, sign(//) is 
considered a complete potential theory of nature with well-defined low energy (weak scale) 
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masses, couplings, etc., obtained from running the masses and couplings from the unification 
scale to the weak scale. All of these low energy parameters must be set for each solution. 
Then, events are generated using 2-to-2 cross section formulae [ 181, H> incorporated into 



Pythia, which also generates initial and final state QCD and QED radiation. Each 2-to-2 
process is derived from a pp collision using the CTEQ2L structure functions. The sparticle 
decays are based on 2-body and 3-body formulae for the decay rates [[H], ^] added 
to Jetset, which performs string fragmentation and hadronization. The final output of 
the event generator is a list of "stable" particles and their 4-momenta. These particles are 
then fed into a model detector which smears momenta based on Gaussian energy resolution 
functions and defines jets based on Et towers with </> x r] segmentation. The detector- 
level variables are used to define the experimental quantities upon which kinematic cuts are 
based. 

Our model detector is based loosely on the CDF detector 22]. The calorimeter is 



segmented A</> x A77 = 0.1 x 0.1 with rj coverage to 77 = 4.2. The hadronic energy resolution 
is chosen so that the jet resolution is roughly 0.7 /VE with a degradation to lA/y/E for large 
77. The electromagnetic energy resolution is 0.2/ yE, and the muon momentum resolution 
is -^7 = v^O-OOOQPT 1 ) 2 + (0.0066) 2 . Electrons, muons, and jets are identified for 77 < 
2.5, though the whole calorimeter is used to define the $t vector. Isolated leptons (/) 
are defined as electrons and muons with Ef tTa = J2i Et - E [ f < 2 GeV within a cone 



R = y Ac/) 2 + At] 2 < 0.4 around the lepton. The sum i is over all leptons, photons, and 
hadrons. Jets j are defined with R = 0.6 and E T > 15 GeV. 

Every signal must go through a process where thousands of supersymmetric solutions 
are simulated each with thousands of events. The final result is a set of observables (cross 
sections, decay rates) which we analyze in the form of scatter plots with each point repre- 
senting a possible supersymmetric theory of nature. 

3 Supersymmetric Signals and Backgrounds 

We consider several potential signals and backgrounds for the Fermilab Tevatron pp 
collider at ^/s = 2 TeV in the following sections organized by the particular detection 
signature. In Section 3.3 we have performed a detailed analysis of the effect of constraints 
on the supersymmetric parameter space with the trilepton signal as our example. 
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3.1 Monolepton and Dilepton Signals and Backgrounds 

We have simulated the monolepton signals from Il^l, Xi Xi production and the dilepton 
signals from Xi xf; X2X?! ^l£*l an d ^R^*r production [23, 24, 25, |2^, 27]. Of course any 
monolepton signal must compete against the huge background from W^{— > £^v) production 
which is ~1 nb. For the supersymmetric signals, the energy spectrum of the lepton from the 
2-body decay of the £ or the (2-body or 3-body) decay of the Xi is generally soft because 
most of the energy is carried by the superpartner decay product Xi- I n addition, the classic 
Jfir signature for SUSY does not easily distinguish signal from background. In the end, we 
could find no set of straightforward cuts which could reliably extract a monolepton SUSY 
signal at the Fermilab Tevatron at any luminosity. 

The dilepton signals appear much more promising, since the large background from 
7*(^ £ + £~) and Z{— > £ + £~) can be significantly reduced with an $t cut and an invariant 
mass cut on opposite sign leptons near the Z mass. However, the background from the 
smaller W (— ► £ u)W^{^ £ T u) production is virtually irreducible at a total leptonic cross 
section of 550 fb (without any cuts). We applied a minimal set of cuts and found the 
dilepton signals from X2X?) ^l£\ and £r£\ are extremely difficult to pull out of the W^W^ 
background of roughly ~140 fb (with our cuts). Only the XiX/i signal is visible in a small 
number of solutions with a chargino mass reach of m-± ~80 GeV, ~110 GeV and ~130 

Xi 

GeV for integrated luminosities of 200 pb^ 1 , 2 fb^ 1 and 25 fb~ . No limits on chargino 
masses could be extracted even at an integrated luminosity of 25 fb . It is possible that 
some of the solutions with detectable trilepton signals described in the next section could 
be confirmed by the dilepton signal [27]. 



3.2 Trilepton Signal and Backgrounds 

The trilepton signal from XiX.® production is probably the most promising signal of 



supersymmetry at a hadron collider [28, p9L iZB, 26, 27]. Although the signal is obtained 



from the chargino decaying to one lepton and the neutralino decaying to two leptons as 
in the previous section, the resulting trilepton signature has two clear advantages over the 
monolepton and dilepton signals. First, the correlations induced between the mass and 
mixing parameters in the CMSSM almost invariably output a lightest chargino which is 
mostly a charged wino (lU^) and a second lightest neutralino which is largely a neutral 
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wino (W 3 ). Thus, the coupling at the W^x^x vertex is near maximal in the CMSSM 
since it is the super symmetrized version of the W^W^T^ 3 Standard Model vertex. Second 
and most important, the trilepton signal has few sizeable backgrounds in stark contrast to 
the monolepton and dilepton signals. In this section we explore all aspects of the trilepton 
signal including a full background analysis and a complete examination of the CMSSM 
parameter space. 

There are six sources of physics^ backgrounds relevant to the trilepton signal: WZ, ZZ, 
it, tb + bt, Z + g and W + g. The WZ, ZZ and it backgrounds have been studied before 
(see e.g. Refs. |2f|). Our results are similar, but more general since we study both the 
traditional backgrounds (first three) and also backgrounds from the leptonic decays of b 
quarks including those from g* — ► bb splitting. The latter can be important especially if one 
considers measuring leptons with low pj- ~ 5 GeV, as we do in this study. For the WZ, ZZ 
and Z+g backgrounds, we properly incorporate the effect of 7* interference in Z production 
and use the full 2 — ► 4 matrix elements inherent in the Pythia Monte Carlo generator. We 
also find it is essential to simulate all possible decay channels of gauge bosons that could 
possibly lead to leptons in the final state {i.e., we include the tau decay channels). In the 
end, we find the WZ and ZZ backgrounds are the most difficult to eliminate, since they 
contain real isolated leptons (distinguishable from leptons from heavy quark decay which 
are generally not isolated from hadrons). 

Below, we enumerate the kinematic cuts applied and the physics motivation behind 
them. 

1. Three isolated leptons (electrons or muons) to reduce QCD backgrounds with pt cuts 
of [pj) ; ; pip) > (10 ; 5 ; 5) GeV for the highest (1), next highest (2), and 
lowest (3) pt leptons. p£P > 10 GeV is necessary for triggering. 

2. \mg+g- — mz\ > 15 GeV for opposite sign, same flavor lepton pair invariant masses to 
reduce Z — > e^e T , [i^y^ backgrounds. 

3. mw > 20 GeV for all lepton pair invariant masses to reduce 7*, b — > £ + X. 

4. Opposite sign £ <yl \ £^ (highest and next highest pt leptons) must not be "back to 
back", \4>2 — 0i I < 2.5 to reduce Z — > t^t t background. 

5. Transverse mass required to be mj < 70 GeV to reduce W backgrounds, where tut 
is constructed from ^M's momentum and the $t vector. 

1 Detectors have additional backgrounds from, for example, 7's that "fake" an electron. The dominant 
background where this can occur is for 7*/Z — ► r + r~ — * e's/fi's with a 7 (and I^t) in the final state. We 
find that at a "fake" rejection rate 1CP 4 our set of cuts effectively eliminates this background. 
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Trilepton Background Results 






a 


after cuts (fb) 




Process 


cut 1 


cuts 1-2 


cuts 1-3 


cuts 1-4 


cuts 1-5 


WZ, Wj* 


22.6 


1.3 


1.0 


0.85 


0.38 


ZZ, Zj*, 7*7* 


5.3 


0.22 


0.16 


0.12 


0.09 


tt 


0.42 


0.33 


0.23 


0.17 


0.06 


Z + g, i*+g 


5.0 


0.64 


0.23 


0.14 


0.14 


Total 


33.3 


2.5 


1.62 


1.28 


0.67 



Table 1: Summary of our background studies of the super symmetric trilepton signal. Cuts 
1-5 are described in the text. The final background estimate using our cuts is 0.67 fb. 

Table [l] shows the reduction of backgrounds as a function of the cuts 1-5. In addition, 
we find no backgrounds from W + g nor tb + hi, which would require two heavy quarks to 
decay to isolated leptons. This reinforces our confidence that we need only consider the 
background from a single heavy quark decay to an isolated lepton. After applying all of 
the cuts listed above, our final trilepton background estimate is 0.67 fb for the Fermilab 
Tevatron pp collider at = 2TeV. 

Using this background estimate, we can calculate the smallest SUSY trilepton cross 
section, folded with the leptonic branching ratios and detection efficiency (denoted a x 
BR x EFF), that is detectable above backgrounds. We find the minimum a x BR x EFF for 
integrated luminosities of 200 pb™ 1 , 2 fb -1 and 25 fb -1 at the Fermilab Tevatron to be 25 
fb, 3.0 fb and 0.82 fb respectively, based on the larger of 5 events or the number of events 
required for a 5<r significance above background. 

In Fig. [l] we have plotted the total supersymmetric a x BR x EFF from X1X2 production 
(including all diagrams) versus the lightest chargino (xi) mass for well over 2000 solutions 
spanning the CMSSM parameter space up to m-± < 500 GeV. Note that a tail of CMSSM 
solutions does exist for m~± > 500 GeV, but we do not consider them here. Each symbol 
represents one solution defined by mo, ni 1 / 2 , Ab tan/3 and sign(^). Each solution, when 
the masses and couplings are run from the unification scale to the weak scale, has its own 
well-defined masses, mixings, branching ratios, etc. The chargino mass represents one of 
these well-defined weak scale observables that is directly related to the production cross 
section (hence, our choice of x-axis). In addition, the branching ratios of the superpartners 
also represent weak scale observables that are crucial to a correct calculation of the trilepton 
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Figure 1: Total super symmetric trilepton signal {a x BRx EFF) after cuts versus the lightest 
chargino mass in the CMSSM. The branching ratio (BR) is defined as the fraction of x^X^ 
events that decay to 3 leptons. The efficiency (EFF) is defined as the fraction of 3 lepton 
events that pass the cuts described in the text. The minimum detectable a x BR x EFF 
for integrated luminosities of 200 pb _1 , 2 fb -1 and 25 fb -1 is shown by the dark horizontal 
lines at 25 fb, 3.0 fb, and 0.82 fb respectively. The different symbols refer to solutions 
showing interesting behavior where the second lightest neutralino (x°) has (A) a neutral 
"invisible" branching ratio (generally X2 ~~ * vjj> then vi, — > x{v) > 90%, (B) a large 
destructive interference in 3-body leptonic decays defined by interference < 0.1 (see Fig. ^), 
(C) a branching ratio to Higgs > 50% dominates, or (D) all other solutions. 
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signal. All decay channels of all possible superpartners that could be produced directly (xf, 
X2), or as a result of decays (gauginos, sleptons and Higgs), must be computed for each 
solution. 

The second lightest neutralino's leptonic branching ratio can sometimes be small enough 
to prevent a detectable trilepton signal even at low neutralino (or chargino) masses. Here 
we elaborate on the neutralino {xq) branching ratios while examining in detail the solutions 
identified in Fig. |l] that have a small leptonic branching ratio. (A) The X2 often decays 
predominantly to neutral ( "invisible" ) products through the 2-body vjj> channel followed by 
vl — > Xi^j or directly through the 3-body Xi u ^- We have identified these solutions in Fig. [I] 
with separate symbols to clearly show they are a source of some low a x BR x EFF solutions. 
(B) The 3-body leptonic decays of the X2 can have destructive interference among the 
diagrams mediated by the £l, (.r and the Z . This interference effect has also been observed 
in Ref. [27]. Once again, we identify such solutions in Fig. [IJwith a separate symbol. Note 
that significant destructive interference occurs only for chargino masses m-± < 200 GeV. 
Generally, the interference is characterized by m| , ra^ ~ 200 GeV. At lower slepton 
masses, it is the slepton mediated 3-body decays that dominate the decay width of X2 
into leptons. At higher slepton masses, the Z mediated 3-body decays dominate. This is 
demonstrated in Fig. where the ratio of the squared-amplitude with and without the 
interference terms for X2 ~~ *■ Xi e+e ~ (which is representative of all the leptonic channels) is 
plotted versus the right selectron mass for all solutions described above. We find destructive 
interference for both signs of fx. The interference for fi > is largest for tan (3 > 10, while 
the interference for fi < occurs over the entire range of tan/3. Finally, (C) the decays of 
X2 Xih° have been explicitly identified in Fig. [I]. We note that to compute the decay 
widths of h°, one must evaluate all masses (and couplings) at the scale of m h o fl30|| . Since 



there are completely general upper mass bounds [31] on m h o <C 2nit, the only decay channels 
open to h° are to bb, t + t~ and lighter particles. When mj, and m T are evaluated at m^o, 
the leptonic decay width of h° — > t + t~ is about 10%. This has a dramatic effect on the 
leptonic branching ratio of the X2 w hen both r's decay leptonically to e's and /x's. Further, 
it is likely that 1-prong pion decays from the Higgs are measurable and could be added to 
the branching fraction; thus all the x's in Fig. |l] (which have a large branching ratio to h°) 
would move up in cross section by a factor of about 2-3. However, accurate estimates of 
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Figure 2: To illustrate the effects of interference in the 3-body decays X2 ~^ Xi^ + ^~^ we 
have plotted the ratio (interference) of the full Z, eR and ei squared-amplitude to the sum 
of the squares of the three individual diagrams for the branching ratio to electrons (which is 
also representative of the branching ratio to muons and taus) against the mass of the right 
selectron. Note that interference is expected to be ~1 if the magnitude of the interference 
is small. 

1-prong backgrounds are more difficult than leptonic backgrounds, so to be conservative 

we do not include isolated pions in our signal. The final result is that we find the trilepton 

signal is present and detectable given a high enough integrated luminosity even for large 

neutralino masses (fn^o ~ m ^ ± <^ 250 GeV) where the x\ ~^ Xi^° decay channel is open. 

Once the cross section and branching ratios are computed, a full event -level simulation 

was performed using the cuts described above. The cuts undoubtedly have an impact on 

the detectable signal, and we quantify this by defining the "detection efficiency" as the ratio 

of the number of 3 lepton events that passed our cuts to the total number of 3 lepton events 

from X1X2 production. We find that the efficiency can vary dramatically over the parameter 

space, as is illustrated in Fig. [| where we plot the efficiency versus the chargino mass. The 

structure of the efficiency plot is clear: At low chargino masses (m~± £ 150 GeV) most of 

the Xi i X2 decays are 3-body with a gradual rise in efficiency due to the increasing energy 

of the leptons. The low efficiency solutions (EFF < 10 -2 , where we found the EFF to 

be > 10 -4 for all solutions) occur when m-± ~ m UL and the 2-body sneutrino modes are 

Xi 

open and there is little energy for the lepton. The neutralino also has analogous problems 
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Figure 3: Detection efficiency of the trilepton signal after cuts versus the lightest chargino 
mass in the CMSSM. The efficiency is defined as the fraction of 3 lepton events that pass 
the cuts described in the text. The different symbols refer to solutions showing interesting 
behavior where the second lightest neutralino (X2) nas (A) a neutral "invisible" branching 
ratio (generally — * ^LM then — ► x\ v ) > 90%, (B) a large destructive interference in 
3-body leptonic decays defined by interference < 0.1 (see Fig. ||), (C) a branching ratio to 
Higgs > 50% dominates, or (D) all other solutions. 
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when to^o ~ mg^ OI ^ R so that 2-body modes are open with one rather soft lepton. At higher 
chargino masses (m~± J> 200 GeV), the 2-body decays Xi ~^ XiW and x\ ~^ Xi% a ^ so 
cause a lower efficiency for the simple reason that our cuts are designed to eliminate real 
W's and Z's. Finally, the somewhat lower efficiency of the x\ — > Xih° mode is visible at 
high chargino masses (m~± > 250 GeV). Note that the efficiencies described here (visible 
in Fig. ^) are automatically included in the a x BR x EFF plot of Fig. [|. 

The results from this trilepton analysis are manifest in Fig. |l]. The Fermilab Tevatron can 
probe chargino masses up to 140 GeV, 210 GeV and 240 GeV with integrated luminosities of 
200 pb" 1 , 2 fb _1 and 25 fb _1 respectively. However, the Tevatron cannot set mass limits on 
charginos or neutralinos from just the trilepton signal. This result follows directly from using 
the full CMSSM parameter space with the complete 2-body and 3-body branching ratios of 
superpartners and a full event -level simulation with realistic detector cuts. The requirement 
for a full simulation of all constrained solutions so that the efficiency is calculated correctly 
for every solution is demonstrated in Fig. ||. Perhaps new experimental or theoretical 
constraints could eventually eliminate the low lying solutions in Fig. |], thereby allowing 
a lower limit on m~± to be set if no signal is found. Nevertheless, we emphasize that a 
significant part of the parameter space can be eliminated if no signal is found. 

3.3 Effects of Constraints on Supersymmetric Parameter Space 

Requiring a constrained parameter space that satisfies our theoretical expectations with 
all current experimental results is very important. To illustrate the effect of constraints 
on the parameter space, we looked at a particular choice of parameters that would give a 
qualitative feel for how b — > 57 and relic density cuts in the CMSSM typically impact the 
trilepton signal. For this example, tan (3 = 5, Aq = 0, and sign(/i) = + (only for illustrative 
purposes in this section), while the values of too and m\i% were selected randomly on a 
logarithmic scale up to 1 TeV. We applied all CMSSM cuts to the data sample except the 
b — > 57 and relic density cuts. 

In Figs. |]a and ^|b we plot the trilepton signal before the b — > 57 and relic density cuts. 
The plotting variables are too (common scalar mass) and mx/2 (common gaugino mass) in 
Fig. |4|a and the more concrete mg fl and to-± in Fig. ^p. (Fig. |]b is nothing more than a 
direct re-mapping of all the points in Fig. |I|a.) The x 's in Figs. |I]a and||b are solutions which 
are detectable through the supersymmetric trilepton signal with an integrated luminosity 



13 



1 1 



ou 



50 100 



500 1000 



m 



1/2 



;GeV) 



50 100 

n>v ( 



;GeV) 



500 1000 



Figure 4: The trilepton signal before b — * sj and relic density cuts. The solutions repre- 
sented by the x's are detectable with an integrated luminosity of 25 fb" 1 at the Fermilab 
Tevatron, while the dots are not. 

of 25 fb^ 1 , while the dots are not detectable without higher integrated luminosity. 

In Figs. [|a and [5|b we produce the same type of scatter plots, this time using the b — > s^y 
observable. All solutions marked with an x are now just those with the additional constraint 
Br(6 -> 57) < 5.4 x 1CT 4 . We choose 5.4 x 1CT 4 as our discriminant value since that is the 



Vo upper limit value published by CLEO |32|. CLEO has also recently reported |S3| a 
measurement of 

Br(6 -> 87) = (2.32 ±0.67) x 10~ 4 . (2) 

This would suggest that we could add 0.67 x 10~ 4 to the central value and declare this our 
la upper bound on Br (6 — ► S7), and then insist that all CMSSM solutions have a predicted 
value of Br(6 — > 57) < 2.99 x 10 -4 . However, the QCD uncertainties [34| in the Br(6 — ► 57) 
calculation (approximately 25%) require us to use a higher "calculation upper limit" than 
the experimental upper limit; therefore, to be conservative we accept all solutions with a 
calculated Br (b — > S7) < 5.4 x 10 -4 . A significant number of solutions in Figs. |5|a and ||b are 
cut out by this constraint. Interestingly, many of the solutions cut out by this constraint 
overlap with light chargino solutions which were not detectable by the trilepton signal. 

The relic density constraint is also an important cut on the CMSSM parameter space. 
We expect that nature is described by an R-parity conserving supersymmetric theory JlC| , 
plf l, and thus the lightest supersymmetric particle (LSP) is absolutely stable. In the early 
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Figure 5: The effect of b — > cut on the parameter space. The solutions represented by 
the x's satisfy the cut Br(6 — > S7) < 5.4 x 10 -4 , while the dots do not. 

universe these stable particles were in thermal equilibrium with the photons until the ex- 
pansion rate of the universe became roughly equal to their annihilation rate. When this 
happens, the LSPs fall out of equilibrium with the photons and their relic abundance stabi- 
lizes. If this decoupling occurs too soon (weak annihilation rates) then the universe becomes 
matter dominated too early. Most observational data indicates that the universe is more 
than about 10 billion years old, and this translates into an upper bound on the relic density 
of LSPs. Quantitatively this upper bound, often called the "age of the universe constraint" , 
is most often expressed clS cL condition on VtLSph ■ 

^Lsph 2 < 10 (age of the universe constraint) (3) 

where h is the Hubble parameter. As part of the CMSSM we require Qj^sph 2 < 1.0. Except 

the age of the universe constraint will exclude 



for a few possible pathological cases [35, 
all SUSY models with large squark and slepton masses |35], [8| . This can be understood by 
realizing that the LSP is mostly (but not completely) bino in the CMSSM, and the bino has 
no coupling to the Z boson and therefore no annihilation channels through an s-channel Z 
are accessible. Then the LSP must annihilate through a t-channel squark or slepton: 



{(JVrel) OC —r. 



m-. 



f 



And since 



LSP 



h 2 oc 



1 



(vv re l) 



then 



flLSph 2 oc r?V~. 



(4) 



(5) 



15 



1000 
500 



> 


100 


CD 


a 


50 


o 




E 


10 



> 

CD 

o 



1000 
500 



lL 



100 

50 




50 100 

m i/2 



500 1000 

;GeV) 



50 100 

m 7+ (GeV) 

A 1 



500 1000 



Figure 6: The effect of the ViLSph? > 0.05 cold dark matter cut and fi^sp/i 2 < 1-0 age of 
the universe constraint cuts on the parameter space. The solutions represented by x 's have 
0.05 < Q, LSP h 2 < 1.0 while the dots do not. 

Therefore, the relic abundance Qlsp grows with the supersymmetric breaking scale. This 
cutoff on supersymmetric masses depends on all the input parameters but is generally 
around 1 TeV. 

The importance of this robust cosmological requirement on studies of supersymmetry 
detectability cannot be overstated. Any successful constraint on the mass of supersymmetric 
particles is obviously of great relevance to a collider program trying to discover or rule out 
supersymmetry. It is for this reason that we have ensured that this constraint is incorporated 
into all solutions analyzed in this paper. 

A lower bound on flLSph 2 can also be obtained from the requirement that the LSPs 
constitute a significant amount of the universe's mass fraction to be a viable cold dark 
matter candidate. Using this requirement, we obtain a quantitative lower bound, 



^LSph 2 <; 0.05 (cold dark matter constraint). 



(6) 



This cold dark matter constraint is not applied to the CMSSM solutions anywhere else in 
the paper except in Figs. || and |?], since no experiment has yet confirmed the identity of the 
cold dark matter. However, we do consider it a major success of the unified supersymmetric 
theories that a stable weakly interacting massive object with a large relic abundance is 
predicted in accordance with astrophysical observations. 

In Figs. 0a and 0b the x's represent solutions which are consistent with both the age 
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Figure 7: The trilepton signal is plotted as in Fig. |] except all solutions (x's and dots) must 
satisfy the b — ► 57 cut and relic density cuts described in the text. Note the large regions of 
parameter space excluded by these cuts, and also the remaining distribution of detectable 
solutions (represented by x's for an integrated luminosity of 25 fb" 1 , as in Fig. 0). 



of the universe and the cold dark matter constraints, summarized by 0.05 < flLSph 2 < 1.0. 
The x's represent all solutions which pass the relic density cut, and the dots represent all 
solutions which lie outside the cut. The dots in the lower left corners of both Figs, ^a and ||b 
are solutions with Q^sph 2 < 0.05, and thus are not interesting dark matter candidates. The 
dots with high mo and high mg R are solutions which have flLSph 2 in qualitative agreement 
with the argument given above that large SUSY scalar masses yield large f^sp/i 2 . Some 
solutions at high mo survive at m^ 2 ~ 100 GeV where the LSPs annihilate through the Z 
resonance. Although the coupling to the Z is quite small, the resonance effect is dominant 
here and the relic density can stay small. 

Finally, Figs. |^a and ^b are identical to Figs. |||a and ||b except now the Br(6 — > 57) 
and £lLSph 2 cuts are included. Many more models are ruled out from the additional cuts, 
and the resulting parameter space provides a more constrained set of solutions that the 
Fermilab Tevatron could either detect or rule out (represented by the x's). The qualitative 
shift in the parameter space by including all constraints on the CMSSM demonstrated in 
this example points to the importance of including all relevant constraints simultaneously 
on supersymmetric solutions for a realistic study of collider capabilities. 
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Figure 8: The correlation between the squark (q) and gluino (g) mass is shown for all 
CMSSM solutions. We have plotted q = ul to be representative of all the squarks. Note 
that relatively few solutions have rriq = nig or nig nig. 

3.4 Missing Et + Jets as a Signal for Squarks and Gluinos 

A classic signature of supersymmetry is multi-jet events with a large $t p7| , 
signal can result from many supersymmetric parton-level processes, and we have simulated 
only a practical subset of these: squark/gluino pair (qq, gg), squark+gluino (qg), and 
squark/gluino + chargino/neutralino production (qx, gxj- In the special case of ~S> niQ 
then nig ~ nig, otherwise nig < niq and we expect the dominant signal to be gg production 
with 3-body decays into other gauginos and jets in that region (see Fig. ||). It is expected 
that the background is dominated by QCD multi-jet production folded with the intrinsic 
]pT resolution of the detector. To estimate the background contribution, we generate all 
QCD parton-level processes for > 30 GeV with initial and final state QCD radiation. 
The IpT resolution is approximated by the CDF-like formula = 0.7\/'SEt, where 
YiEt is the scalar sum of all the jet transverse energy. The resolution degrades with 
increased jet activity, so only high pt events can generate a large ^t- The contribution of 
each event is folded with the probability that the measured $t fluctuates to $^ ut or more 
using the Gaussian formula: 

P(E T > $£ ut ) = -= / e _t dt where t' = (7) 

V 7T Jf OA y/2 



38, 



This 
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The resultant weighted cross section for E^ ut = (0,50,75,100,125) GeV is (0.13 fib, 1.3 pb, 
2 fb, ~ 0, ~ 0). However, a Gaussian approximation to the resolution is not realistic, since 
non-Gaussian tails are known to be important. To test the sensitivity of our estimate, we 
have added a term to our probability distribution so that the total probability of mismea- 
surement from 2.58-6.00cr is still ~ 1% but the probability is evenly distributed. For this 
case, our background estimate is (0.13 fib, 115 pb, 1.8 pb, 103 fb, ~ 0). Because of the 
strong sensitivity to the non-Gaussian tail, we choose $^ ut = 75 GeV. This reduction of 
the background is conservative since we have not utilized the altered kinematics resulting 
from the loss of jet energy. We further apply a transverse sphericity cut St > 0.2 to reduce 
this pure jet background to a negligible level. The other physics backgrounds we considered 
all have $t from Z{— > vV) + g, Z(— ► r + r~), W{-^ Iv) + g, and tt. Since we veto events 
containing isolated muons and electrons with px > 15 GeV, W(— > rv) is a large potential 
background. Backgrounds are also reduced by requiring that the sum E^ + Ej? > 300 
GeV, where j\ and j'2 are the two highest Ej- jets. Finally, we require that A(p > 0.5 
between each jet and the $t direction to reduce the fake $t from energy lost in cracks or 
shower fluctuations. Using these cuts, the cross sections are 5 fb, 11 fb, and 24 fb for the 
Z, W, and tt backgrounds respectively. 

After a full simulation we find the gluino mass reach for integrated luminosities of 200 
pb" 1 , 2 fb" 1 and 25 fb" 1 to be ~300 GeV, ~350 GeV, and ~400 GeV respectively, at 10a 
signal significance (see Fig. ||). Here significance is defined as the number of signal events 
divided by the square root of the number of background events, for events that pass the 
above cuts with = 75 GeV. Of course, this is not a useful procedure unless one has 
reasonable knowledge of the normalization of the experimental background I^t distribution 
as well as all the real backgrounds, which we assume will occur. Note that we do not make 
explicit use of the signal or background shape, just the number of events above ^^ ut . 

We have studied several kinematic variables to see if the resulting distributions can 
make the signal more convincing or help us extract any information about the superpartner 
masses. As an example, we choose a model with rrig = 298 GeV and rriq = 311 GeV that 
has a total production cross section of 1.4 pb. The invariant mass found by summing the 
four vectors of the 4 highest Et jets with the four vector {^t-,^t) is shown in Fig. |n] for 
the signal (dark upper bins) and the backgrounds (lower light or hashed bins). The inset 
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Figure 9: The a x EFF is plotted for the I^t + jets signal versus the gluino mass. The 
efficiency (EFF) is defined as the fraction of events that pass the cuts described in the 
text. Horizontal lines indicate the gluino mass reach of the Fermilab Tevatron at integrated 
luminosities of 200 pb -1 , 2 fb -1 , and 25 fb -1 . The dotted line represents the background 
cross section with our cuts. The vertical banding in gluino mass is due to numerical sampling 
and is not physically significant. 



20 



70 - 




M m (GeW) 



Figure 10: The invariant mass distribution (Mjjjj x ) found by summing the four vectors 
of the 4 highest Et jets with the four vector ($t,$t) is plotted for an example model 
with rrig = 298 GeV and rriq = 311 GeV. The dark shaded region on the top is the signal, 
while the lower three shaded regions are the backgrounds from Z, W and it production 
respectively. The inset figure shows the corresponding $t distribution for the same model. 

figure shows the distribution for the same processes. In the invariant mass spectrum 
the signal clearly peaks around 2rrig ~ 2m,q ~ 600 GeV, while there is no such peak in the 
]fir distribution. Since the backgrounds from Z, W, and it peak at lower invariant masses, 
one can see how the invariant mass distribution provides an effective means to separate the 
signal from background. 

3.5 Signals of Stop Production 

The signal from stop production has been known for some time to be quite promising [|3S| , 
p3| , |40"| , p| |42"[ since the lightest stop mass eigenstate (t±) is generally lighter than all the 
other squarks. This follows from the substantial mixing between the top squark weak 
eigenstates ti, and tn caused by a large mixing term in the stop mass matrix proportional 
to the top mass. Furthermore, the large top Yukawa coupling in the running of the stop 
scalar mass terms reduces the mass of the stop with respect to the other squark masses. 
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The result is a large mass splitting between the stop mass eigenstates t\ and ti- Hence, t\t\ 
is more readily pair produced than the other squarks. The decay i\ — > x[b dominates if it 
is kinematically allowed, yielding a final state analogous to t — ► W + b in the SM. For those 
models of the CMSSM where this decay is not allowed, we assume the decay t\ — > x? c with 
a 100% branching ratio. 

A stop search is similar to a top search if t\ — > though we take advantage of the 
different kinematics to separate stop from the top background. The full set of backgrounds 
considered is it, tv) + g, + 5, W ± W T , and t + t~). Additional jet 

activity is generated by gluon splitting and initial and final state radiation. We classify 
the events into 3 channels: 1) dilepton, 2) W + jets, and 3) dijet. Channels 1) and 2) are 
distinct from a top signal because the lepton (W + jets) or leptons (dilepton) have smaller 
Pt and an Et inconsistent with W decay. Dijet events 3) are Et + jets events with no 
isolated leptons. We identify dilepton events by the following kinematic cuts: 

1. Two isolated electrons or muons with Et > 10 GeV. 

2. Two or more jets, one or two 6-tagged. 

3. E T > 25 GeV; if 25 < $ T < 50 GeV then an additional cut A0 > 0.3 is applied, 
where Acf> is measured between $t and any of the jets. 

4. \m u > — mz\ > 15 GeV, and m u i > 10 GeV for leptons £, ( . 

5. m^ AX < 100 GeV, where m£ I fe AX is the highest invariant mass between the leptons 
and the 6-jets. 

W + jets events are defined by: 

1. One isolated electron or muon with Et > 10 GeV. 

2. Two or more jets, one or two 6-tagged. 

3. S T > 25 GeV. 

4. rriT < 75 GeV, where the transverse mass rriT is constructed from the lepton momen- 
tum and the $t- 

5. mf b AX < 100 GeV, where 

m £fe AX i s highest invariant mass between the leptons 

and the 6-jets. 
Finally, dijet events are classified by: 

1. Two or three jets with E T > 20 GeV, and no isolated leptons. 

2. E T > 75 GeV. 
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3. A(j) jlj2 < 3. 

4. EE 3 T < 150 GeV. 

5. Tnj lt j 2 > 120 GeV, where mj 1 j 2 i s the invariant mass of the highest (ji) and next 
highest (j 2 ) E T jets. 

For event types 1) and 2), where a b- tag is required, we assume an efficiency ej, per single 
tag independent of pt or rj of the tracks. Using e& = 0.3, the probability to tag one or two 
6's is then roughly 0.5. For dijet events, it was shown previously that soft lepton tagging of 
the c-jet does not dramatically improve the stop search EJ, so we do not consider it here. 



Our results are illustrated in Fig. |11|, where the cross section folded with the detection 
efficiency (a x EFF) is plotted versus the stop mass for the dilepton, W + jets, and dijet 
channels. We find that for an integrated luminosity of 2 fb _1 , the stop mass reach is up to 
~160 GeV at 5a significance in the W + jets channel. For 25 fb , the stop mass reach is 
up to ~200 GeV at 5a significance in the W + jets channel and the dilepton channel. The 
few models with i\ — ► \\c are visible in the upper left corner of the dijet channel graph. 

3.6 Comment on Top Decays to Stop Signals 

Another interesting signal of top squarks is through the decay of the top quark into a 
top squark and an LSP (t — ► tiXi) [39, We have analyzed all of the CMSSM solutions 



to find those which kinematically allow this decay, and which also change the branching 
fraction of t — * W + b significantly (however, see Section 4). All solutions which satisfy the 
kinematic requirement also give Br(£ — ► tiXi) ~ 15%. Since the x? is invisible and the 
stop will be hard to see with the soft jets originating through t\ — > x± c decays, perhaps the 
only way to detect the effects of the light stops and LSPs is through modification of the 
Br(i — ► W + b) ~ 1 branching ratio. 

The best way to test for the existence of these stops is through stop production. Since 
in the mt > + m^,o region of parameter space the stops will decay primarily through a 
one-loop diagram to a charm and an LSP, good charm tagging would be extremely helpful 
if this signal has a chance of being distinguished from other copious jet sources. 

4 Light Charginos and Stops Beyond the CMSSM 

All analyses in this paper are based on the constrained minimal supersymmetric Stan- 
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Figure 11: The a x EFF for stop production is plotted versus the stop mass in the dilepton, 
W+jets and dijet channels. The efficiency (EFF) is defined as the fraction of events that 
pass the cuts described in the text for each channel. Horizontal lines indicate the reach of 
the Fermilab Tevatron at integrated luminosities of 200 pb _1 , 2 fb _1 , and 25 fb -1 . The 
dotted lines represent the total background cross section for each channel with our cuts. 
Note the total stop production cross section (with no cuts) is also shown for comparison. 
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dard Model (CMSSM) described in the introduction. The CMSSM is more general than the 
SM since it includes the SM fully, it is a consistent theory, it incorporates phenomena not 
explained by the SM such as the apparent unification of the gauge couplings, and provides 
a derivation of the Higgs mechanism rather than assuming it as in the SM. Further, we 
have seen from Section 3 that the CMSSM makes a variety of testable predictions. 

However, there is a hint that the CMSSM may not be entirely correct because its the- 
oretical assumptions about super symmetric parameters are too restrictive. The purpose of 
this section is to point out this possibility and the implications for supersymmetry searches 
at the Fermilab Tevatron. 

There are two effects that are clues to go beyond the CMSSM. First, the Br(Z — ► bb) is 
larger than its SM value by about 2a. Second, the value of a s deduced in LEP analyses from 
the Z width is larger by about 2a than the value of a s deduced other ways. While neither of 
these effects is of great "statistical" significance, both are based on a number of independent 
measurements over several years from several detectors; the errors are mostly systematic 



and theoretical. The possible importance of these effects is stated in Ref. [17], which finds 
that by including light superpartners in the analysis of the LEP data, both effects can be 
simultaneously explained. In particular, it is remarkable that the LEP a s measurement 
becomes consistent with other a s measurements (at a s ~ 0.112) when light charginos and 
stops are included. Further, Ref. [|17| reports a global analysis of the LEP/SLC data to 
ensure that no other observable is adversely affected by the light superpartner contribution 
(such as Ap or raw). One finds the supersymmetry fit is actually better than the SM 
fit, giving better agreement to, for example, the Ri, and Alr measurements. The better 
agreement requires both m~± and rar to be less than about 100 GeV. 

Surprisingly, the properties of the chargino and stop are sufficiently well determined by 
the above analysis that one is forced into concluding that the theoretical assumptions of the 



CMSSM cannot be fully correct. The required mixture [15| of gaugino and Higgsino in the 
chargino, and the left-right mixing for the stop, are not allowed by the CMSSM, though 
they are fully allowed by the theory when particular assumptions about parameters in the 
CMSSM are relaxed. There has not been time to construct constrained models based on 
the non-minimal theory, but we can summarize the expected impact on the opportunity to 
detect superpartners at the Fermilab Tevatron. 
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First, and most important, the light chargino and stop should be fully produced at the 
Fermilab Tevatron. The stop production is predominantly via gluons and is unaffected since 
the cross section is independent of II, Ir mixing. The chargino production rate will also 
not be greatly affected over the parameter space. Kinematically, since charginos and stops 
are lighter than about 100 GeV, there is no suppression at the Fermilab Tevatron. 

The signatures and detection will be different, and can be more difficult. The analysis 
for stop production and detection is largely the same. However, charginos that resolve the 
a s "crisis" are more Higgsino-like than CMSSM charginos, which implies the LSP mass 
is closer to the chargino mass, and so the resulting leptons are softer (hence fewer leptons 
pass the pr cut). One new helpful feature does enter, in that now all stops are lighter than 
tops, so detecting stops in top decay is a major opportunity. The Br(t — > t + LSP) has 
a factor of Gi?mf and the usual phase space factor, plus a factor of iV^/sin 2 f3 that is the 
supersymmetric modification from the Higgsino wave functions and the top-stop coupling. 
Numerically, we find 0.15 < Br(t -► t + LSP) < 0.6 for m-± ~ 60 GeV, to 0.05 < 
Br(t -> t + LSP) <> 0.3 for m-± ~ 90 GeV. 

Once the decay t — > t + LSP has occurred, the analysis of t decays is as before, with 
t\ — > xfb dominating if it is open, and t\ — ► Xi c otherwise. For chargino detection a new 
analysis is required to correctly calculate the BR x EFF for the trilepton signal, and the 
needed models have not yet been analyzed. Qualitatively, as remarked above, the leptons 
will be a little softer, so lowering the lepton triggers will increase the efficiency, and gains 
by going to softer leptons may be major. The interference effect and the invisible vjjv mode 
may not cause such a large suppression of the leptonic branching ratio of the neutralino as 
in the CMSSM, possibly improving the situation for some of the parameter space. 

In summary, if the LEP R\ } measurement and the a s crisis are telling us that light 
charginos and stops exist, then they are already being produced at the Fermilab Tevatron. 
The challenge is to detect them. The analysis based on the CMSSM, or an extended theory, 
can help thinking about signature and detection efficiency, and help the interpretation if no 
signal is found. If light charginos and stops are there, high luminosity and a careful search 
is needed to find them. 
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5 Conclusions and Comments 



We have conducted a thorough study of many possible supersymmetry signals at a 
luminosity upgraded Fermilab Tevatron collider. Our studies have complemented and of- 
ten improved on others using our full event-level Monte Carlo based on Pythia/Jetset 
coupled with a comprehensive approach to constrained minimal supersymmetric parame- 
ter space. We conclude, consistent with earlier work, that the trilepton signal and gluino 
production are the most useful signals to discover supersymmetry at the Fermilab Tevatron 
collider. Although no limits could be placed on the mass of the lightest chargino if no sig- 
nal were detected, the discovery potential at Fermilab is quite large, extending to chargino 
masses well beyond that which can be probed by LEP II. With integrated luminosities of 
200 pb _1 , 2 fbr 1 and 25 fb _1 , the Fermilab Tevatron can reach to chargino masses of 140 
GeV, 210 GeV and 240 GeV respectively. Gluino searches at Fermilab will be able to detect 
solutions with m~ g up to ~300 GeV, ~350 GeV, and ~400 GeV to 10cr significance with 
integrated luminosities of 200 pb _1 , 2 fb _1 , and 25 fb _1 . The stop production signal at 
Fermilab has a stop mass reach up to ~160 GeV at 5cr significance in the W + jets channel 
with an integrated luminosity of 2 fb _1 . At an integrated luminosity of 25 fb _1 , Fermilab 
has a stop mass reach up to ~200 GeV at 5a significance for stop production in both the 
W + jets and the dilepton channels. 

While searches at LEP can likely find or exclude chargino and stop masses up to nearly 
the beam energy searches at the Fermilab Tevatron have a reach ~2-3 times that of 
LEP II. But, the Fermilab Tevatron cannot set mass limits if no signal is found because sets 
of parameters exist that give smaller a x BR x EFF than can be detected. However, in the 
absence of a discovery what matters is establishing tighter experimental constraints on the 
parameters; in Section 3.3 we showed how a variety of information will combine to reduce 
the allowed parameter space, so that increasingly unique and testable predictions can be 
made. In this pursuit, LEP II and Fermilab are comparably powerful and somewhat com- 
plementary. Searches for supersymmetry utilizing existing and upgraded collider facilities 
have a great potential for discovery, and if no discovery occurs at a given energy or lumi- 
nosity then the results help sharpen both our understanding of the theory and predictions 
for a variety of other experiments. 
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